Abstract
Recently, growth mechanism of firms in complex business networks became new targets of scientific study owing to increasing availability of high quality business firms’ data. Here, we paid attention to comprehensive data of M&A events for 40 years and derived empirical laws by applying methods and concepts of aggregation dynamics of aerosol physics. It is found that the probability of merger between bigger firms is bigger than that between smaller ones, and such tendency is enhancing year by year. We introduced a numerical model simulating the whole ecosystem of firms and showed that the system is already in an unstable monopoly state in which growth of middle sized firms are suppressed.
Similar content being viewed by others
Introduction
The analysis and modeling of our society is providing profound new insights such as cascading failures1 and resilience2, 3 using the lens of the science of complex systems4,5,6,7. In earlier research on a national network of transactions between businesses, which can be regarded as a typical real world example of a complex network having both the basic properties of small world8 and scale-free9 accompanied with various scaling relations10,11,12,13,14,15,16,17,18, it was found that the balance between new entry, bankruptcy and merger plays the key role in the time evolution of this system19, 20. Here, we focus on the recent comprehensive data of Mergers and Acquisitions (“M&A”) in Japan, more than 41 thousands events, and find that the network evolution is essentially unstable due to two basic inherent business dynamics; the preferential selection, whereby acquiring companies tend to be bigger than their targets, and concentrative M&A, in which the probability of merger among bigger companies is higher. As a result, the time evolution equation with the parameters estimated from recent data indicates that there is no stable steady solution in the size distribution of firms. Our numerical simulation shows that the selective M&A mechanism leads to significant concentration resulting in the appearance of business monopolies and oligopolies, akin to those types described in other research such as super connected companies21. Importantly, such companies suppress growth of middle size firms, and therefore suppress the emergence of new challengers that are fundamental to the functioning of an ecosystem. Sooner or later when the giant fails, a turbulent period follows until another giant appears repeating with a typical time scale of a few tens of years.
Results
Increasing of M&A activities
According to comprehensive business firms data (see Supplementary Information; Data source), the annual number of M&As in Japan in the last 30 years has increased more than 10 times as in other countries such as the EU and China in Fig. 1a (all data of the whole world, the US, the EU and China are sourced from the Institute for Mergers, Acquisitions and Alliances). This tendency is clearly illustrated by an actual merger history diagram of a particular firm, which is the center position node (Fig. 1b). M&A activity has a potent influence on business transactions and hence is one of the most important corporate activities. We next focus on inter-firm transactions network aspect of M&A. Figure 1c shows the spatial distribution of firms sized by number of transaction partners, or link number, k within Kyoto Prefecture in 1994 and 2014. In this case, it appears that majority of firms with high link number in 2014 have been involved in mergers (yellow circles), although it is equally noticeable that the majority of firms in 1994 had no merged ancestors (white circles).
Scaling behaviour and power laws
It has been shown mathematically that M&A activity is one of the main causes of the power law link number distribution of firms with exponent 1.4 measured by k 19. Firm sizes can equally well be measured by number of employees E, by annual sales S, and number of ancestors A correlates with number of employees and annual sales. There are scaling laws between pairs of k, E, S, and A such that \(S\propto {k}^{1.3}\), \(S\propto {E}^{1.3}\), \(E\propto {k}^{1.0}\), \(k\propto {A}^{1.4}\) (Supplementary Fig. 2), and we apply k as the measure of firm size in the following discussion (it can be substituted E, S or A given that there are similar scaling laws among those quantities as previously described). Moreover, all these quantities follow power law cumulative distributions with exponents 1.4, 1.0 and 2.4, respectively (Supplementary Fig. 1).
Statistical characteristics of M&A activities
By making use of the data, we are able to determine the statistical characteristics of firms M&A activities. Firstly, we compare a merged entity with the sum of an acquirer firm a and target firms t (Fig. 2a). As a result, we find that M&A almost keeps the conservation law for each volume such as k, E and S (Fig. 2b and Supplementary Fig. 3). Based on the simple law, we can draw inspiration from Smoluchowski coagulation equation22 to observe a frequency of merger \(K({k}_{a},{k}_{t}){c}_{{k}_{a}}(\tau ){c}_{{k}_{t}}(\tau )\) and merger kernel K(k a , k t ). Note that k a is a link number k of an acquirer firm a, and k t is that of a target firm t, and c k (τ) is a probability density of firms with k in time τ. We then observed \(K({k}_{a},{k}_{t}){c}_{{k}_{a}}(\tau ){c}_{{k}_{t}}(\tau )\) (Fig. 2c) and K(k a , k t ) (Fig. 2d). With regards to the process of mergers and acquisitions among firms, we identify two key characteristics;
-
1.
Acquirers tend to be bigger than their targets and to have a preference (Fig. 2c).
-
2.
The probability of merger between bigger firms is higher than that between smaller ones (Fig. 2d).
We simultaneously find the K(k a , k t )’s scaling laws with k a and k t (Fig. 2e,f) as follows;
Here, \((\alpha \mathrm{,\ }\beta )\simeq (0.8\pm 0.1,\mathrm{\ 0.6}\pm 0.2)\) from 1993 to 1995, and \((\alpha \mathrm{,\ }\beta )\simeq (1.1\pm 0.2,\mathrm{\ 0.7}\pm 0.3)\) from 2013 to 2015. We make a distinction between two separate elements of the M&A process as \({k}_{a}^{\alpha }\) corresponds to the willingness of a company to merge, whereas \({k}_{t}^{\beta }\) is the choice, or the selection, of a given partner. We observed these exponents by the least squares method with coefficient of determination R 2 ≥ 0.9 and the error-bars show standard deviation. That is, these exponents characterize the tendencies toward M&A activity, which has statistically significantly increased recently. We note that we have accumulated 3 years of data (1993–1995; 2013–2015) in order to accurately estimate the statistical properties of the merger kernels of 1994 and 2014 by increasing the sample sizes. Besides, as for the established firm case in 1994 and 2014, we can observe certain level of preferential attachment9, 19. Supplementary Fig. 4 shows a distribution \(\kappa (k)=\frac{Q(k)}{N(k)}\sim {k}^{\lambda }\) where Q(k) is the probability of a new entrant connecting to an old firm of size k, N(k) is the number of firms with k, and λ is the exponent of preferential attachment23. Each exponent in 1994 and 2014 is about \(\lambda \simeq 1.0\pm 0.1\) estimated by the least squares method with coefficient of determination R 2 ≥ 0.9 and error-bars show standard deviation, which is statistically consistent.
Numerical analysis of selective M&A mechanism effects on firm-ecosystem
In order to better understand, we introduce simulation based on network evolution model and numerically analyze our ecosystem’s time evolution by using inter-firm business transactions networks. There are various types of models that focus on network evolution9, 24,25,26,27. Among them, Miura et al.19 introduced so-called MTT model and found that coagulation, namely M&A, is one of the important processes to reproduce distribution of transaction partners k that follows a power law consistent with that of Japanese inter-firm business transactions network (Supplementary Fig. 1). Furthermore, their model considers firm events such as new establishments and bankruptcies, and we decided to generalize MTT model to replace the merger kernel by the empirical form \(K({k}_{a},{k}_{t})\propto {k}_{a}^{\alpha }{k}_{t}^{\beta }\) in equation (1).
Based on this revised model, we conducted Monte Carlo simulation to analyze time evolution of the ecosystem of business firms numerically (see Methods and Supplementary Fig. 5). We start with N 0 = 100,000 firms without any business transactions and evolve the system by choosing one of the three events, new establishments, M&A and bankrupts, stochastically for 10,000,000 time steps where 8,000 steps corresponds to one year. Here, occurrence probabilities of new establishments, M&A, and bankrupts are 0.5, 0.1, and 0.4, respectively, estimated by the state of real ecosystem in Japan. The parameter-sets (λ, α, β) are empirically observed (1.0, 0.8, 0.6) and (1.0, 1.1, 0.7) for the year of 1994 and 2014, respectively (Fig. 2e,f). The difference of coagulation kernel causes a big difference in stability of systems as explained below.
We observed examples of simulated time series after 2,000,000 steps; 1st, 2nd, 100th (top 0.1%), 1,000th (top 1%), 10,000th (top 10%) and median (top 50%) to check a system’s stability. As for the year 1994 (Fig. 3a), the state of the system is almost stable. By contrast, the state is unstable in 2014 (Fig. 3b); there are larger fluctuations in 2014 than in 1994 because of the occurrences of dominant entities for every few decades. The number of transaction partners of the dominant company grows exponentially with the age of the firm on average (Supplementary Fig. 6c,i). Especially, focusing on a growth speed of nodes which have the most number, top degree nodes have a lead of 10 times over second-top degree nodes by 45 ± 22 years on average after disappearing ex-top degree nodes (Supplementary Fig. 7).
We hence numerically confirm that the selective M&A mechanism that prefers to merge with larger size of firms leads to significant concentration and the dominance of large entities, and a consequent suppress of growth of middle-sized firms (see Supplementary Information; A consequence suppress of growth of middle-sized firms). An example of abrupt changes in the number of transaction partners k of the dominant company is shown in Fig. 3c, which represents time evolution of real firms up to 2014. The firm that had second largest number of transaction partners (green marks) from 2000 to 2005 merged its 27 subsidiaries in 2006 and since then, it has kept the 1st seat (red marks). Figure 3c also shows a speculation of firm-ecosystem from 2015 to 2035 by using parameter-set in 2014, which the initial condition is the real network in 2014 with N 0 = 1,132,629; 1st, 2nd, 1,000th (top 0.1%), 10,000th (top 1%), 100,000th (top 10%) and median (top 50%). Only two superordinates will continue growing.
Discussion
We empirically find that our Japanese business firms have the specific property of the selective M&A mechanism that benefits large-sized firms growth. This means that there is an arbitrary “selection and concentration” process within our ecosystem, which leads to an imbalance to the latter. Within the real world Japanese data, we are also able to identify that the two separate dynamics present within M&A, namely the willingness of a company to merge, and its choice of partner, with the former having a higher influence than the latter. The increasing strength from 1994 to 2014 of the Cumulative Advantage mechanism that underpins the willingness to merge resulted in an acceleration of the inequality of income among companies as observed by the increase in the GINI index28 (Supplementary Fig. 8). Indeed the concentration of income within large companies is documented by other recent researches21, but no mechanism was explained. Therefore, importance must be given to the fact that strong selective and concentrative merger allows exclusive possession of the wealth to a mega-firm deriving the unstable fragile economy by its failure with potentially major negative consequences to the global economic network in the long timespan. We believe that our methods based on big data analysis and modeling can be added to basic scientific tools for stability analysis of complex ecosystems in the real world not limited to the business firms’ ecosystem.
Methods
Monte Carlo Simulation
As for our numerical analysis to estimate our ecosystem’s future behaviour, we conducted Monte Carlo simulation as follows;
-
Start with N 0 firms without any business transactions (N 0 = 100,000 in Fig. 3a,b, and N 0 = 1,132,629 in Fig. 3c).
-
Choose one of the following three events stochastically. The occurrence probabilities of new establishments, M&A, and bankrupts are denoted by r n , r m , and r b , respectively.
New establishments
A new firm having four transaction partners is added; it is roughly consistent with the rate between number of firms and number of transaction partners. Besides, each transaction partner is connected to a firm chosen randomly following the preferential attachment rule with exponent λ (Supplementary Fig. 4).
M&A
Firms are randomly chosen following the merger kernel \(K({k}_{a},{k}_{t})\propto {k}_{a}^{\alpha }{k}_{t}^{\beta }\) to choose an acquirer firm a and a target firm t (Fig. 2e,f). All the transaction partners connected to the target firm t are also rewired to the acquirer firm a.
Bankrupts
A randomly chosen firm is removed, along with all transaction partners connected to this firm. This is because age of firms follow exponential distribution; it is roughly consistent with the simple assumption that a firm disappears randomly following a Poisson process19, 29 (Supplementary Fig. 6a).
-
Repeat Step2 for 10,000,000 times.
We also illustrate our algorithm by schematic diagram in Supplementary Fig. 5b. The occurrence probabilities of new establishments, M&A, and bankrupts are denoted by r n , r m , and r b , respectively, satisfying \({r}_{n}+{r}_{m}+{r}_{b}=1\) (systems are quasi-stationally states when \({r}_{n}={r}_{m}+{r}_{b}\) in the number of firms). We here used \({r}_{n}=0.5\), \({r}_{m}=0.1\) and \({r}_{b}=0.4\) corresponding to the state of ecosystem in Japan (Supplementary Fig. 5c). In the case N 0 = 100,000, 8,000 steps corresponds to one year which is estimated by matching the age of fimrs’ distribution (Supplementary Fig. 6a,g). In the case of simulation of Fig. 3c after 2015, the real network structure is used as the initial condition and the parameters are given by the real values measured in 2014. One year is given by 80,000 steps.
References
Buldyrev, S. V., Parshani, R., Paul, G., Stanley, H. E. & Havlin, S. Catastrophic cascade of failures in interdependent networks. Nature 464, 1025–1028 (2010).
Cohen, R., Erez, K., Ben-Avraham, D. & Havlin, S. Resilience of the Internet to random breakdowns. Phys. Rev. Lett. 85, 4626–4628 (2000).
Gao, J., Barzel, B. & Barabási, A.-L. Universal resilience patterns in complex networks. Nature 530, 307–312 (2016).
Nepusz, T. & Vicsek, T. Controlling edge dynamics in complex networks. Nature Phys. 8, 568–573 (2012).
Barzel, B. & Barabási, A.-L. Universality in network dynamics. Nature Phys. 9, 673–681 (2013).
Cornelius, S. P., Kath, W. L. & Motter, A. E. Realistic control of network dynamics. Nature Commun. 4 (2013).
Majdandzic, A. et al. Spontaneous recovery in dynamical networks. Nature Phys. 10, 34–38 (2014).
Watts, D. J. & Strogatz, S. H. Collective dynamics of ‘small-world’ networks. Nature 393, 440–442 (1998).
Barabási, A.-L. & Albert, R. Emergence of scaling in random networks. Science 286, 509–512 (1999).
Stanley, M. H. R. et al. Scaling behaviour in the growth of companies. Nature 379, 804–806 (1996).
Amaral, L. A. N. et al. Scaling behavior in economics: I. Empirical results for company growth. J. Phys. I France 7, 621–633 (1997).
Takayasu, H. & Okuyama, K. Country dependence on company size distributions and a numerical model based on competition and cooperation. Fractals 6, 67–79 (1998).
Axtell, R. L. Zipf distribution of US firm sizes. Science 293, 1818–1820 (2001).
Takayasu, M. et al. Massive Economics Data Analysis by Econophysics Methods-The case of companies’ network structure. Earth Simulator Center Report 263, 263–268, (2008).
Fu, D. et al. The growth of business firms: Theoretical framework and empirical evidence. Proc. Natl. Acad. Sci. USA 102, 18801–18806 (2005).
Watanabe, H., Takayasu, H. & Takayasu, M. Relations between allometric scalings and fluctuations in complex systems: The case of Japanese firms. Phys. A 392, 741–756 (2013).
Takayasu, M., Watanabe, H. & Takayasu, H. Generalised central limit theorems for growth rate distribution of complex systems. J. Stat. Phys. 155, 47–71 (2014).
Daepp, M. I., Hamilton, M. J., West, G. B. & Bettencourt, L. M. The mortality of companies. J. R. Soc. Interface 464, 20150120 (2015).
Miura, W., Takayasu, H. & Takayasu, M. Effect of coagulation of nodes in an evolving complex network. Phys. Rev. Lett. 108, 168701 (2012).
Viegas, E., Cockburn, S. P. & Jensen, H. J. The dynamics of mergers and acquisitions: ancestry as the seminal determinant. Proc. R. Soc. A 470, 2171 (2014).
Vitali, S., Glattfelder, J. B. & Battiston, S. The network of global corporate control. PloS One 6, e25995 (2011).
Smoluchowski, M. V. & Vortrage, D. Uber Diffusion, Brownsche Bewegung und Koagulation von Kolloidteilchen. Phys. Zeits. 17, 557–585 (1916).
Krapivsky, P. L., Redner, S. & Leyvraz, F. Connectivity of growing random networks. Phys. Rev. Lett. 85, 4629 (2000).
Albert, R. & Barabási, A.-L. Topology of evolving networks: local events and universality. Phys. Rev. Lett. 85, 5234–5237 (2000).
Dorogovtsev, S. N., Mendes, J. F. F. & Samukhin, A. N. Structure of growing networks with preferential linking. Phys. Rev. Lett. 85, 4633 (2000).
Krapivsky, P. L. & Redner, S. Organization of growing random networks. Phys. Rev. E 63, 066123 (2001).
Moore, C., Gourab, G. & Newman, M. E. Exact solutions for models of evolving networks with addition and deletion of nodes. Phys. Rev. E 74, 036121 (2006).
Gini, C. Variabilità e mutuabilità. contributo allo studio delle distribuzioni e delle relazioni statistiche. Tipogr. di Cupini (1912).
Goto, H., Takayasu, H. & Takayasu, M. Empirical Analysis of Firm-Dynamics on Japanese Interfirm trade Network. Social Modeling and Simulation plus Econophysics Colluquium 2014, Proc. Int. Conf., 195–204 (2014).
MATLAB 8.5, The MathWorks, Inc., Natick, Massachusetts, United States, https://www.mathworks.com/.
Arc GIS 10.2, Environmental Systems Research Institute, Inc. (ESRI), Redlands, California, United States, http://www.esri.com/arcgis/about-arcgis.
Acknowledgements
The authors appreciate Teikoku Databank, Ltd., Center for TDB Advanced Data Analysis and Modeling for providing both the data and financial support. This study is partially supported by the Grant-in-Aid for Scientific Research (B), Grant Number 26310207 and JST, Strategic International Collaborative Research Program (SICORP) on the topic of “ICT for a Resilient Society” by Japan and Israel. This research was partially supported by MEXT as “Exploratory Challenges on Post-K computer (Study on multilayered multiscale spacetime simulations for social and economical phenomena)”.
Author information
Authors and Affiliations
Contributions
H.G. performed most experiments, interpreted the results, and co-wrote the manuscript. E.V. performed experiments, interpreted the results, and co-wrote the manuscript. H.J.J. supervised experimental designs, interpreted the results and co-wrote the manuscript. H.T. designed experiments, supervised experimental designs, interpreted results, and co-wrote the manuscript. M.T. conceived the project, designed experiments, supervised experimental designs, interpreted results, and co-wrote the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Goto, H., Viegas, E., Jensen, H.J. et al. Appearance of Unstable Monopoly State Caused by Selective and Concentrative Mergers in Business Networks. Sci Rep 7, 5064 (2017). https://doi.org/10.1038/s41598-017-05362-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-017-05362-5
This article is cited by
-
Assembling real networks from synthetic and unstructured subsets: the corporate reporting case
Scientific Reports (2019)
-
Modeling and simulation of Japanese inter-firm network
Artificial Life and Robotics (2019)
-
Smoluchowski Equation for Networks: Merger Induced Intermittent Giant Node Formation and Degree Gap
Journal of Statistical Physics (2018)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.